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Abstract. We introduce MADE3 (Modal Aerosol Dynam- 
ics model for Europe, adapted for global applications, 3rd 
generation; version: MADE3v2.0b), an aerosol dynamics 
submodel for application within the MESSy framework 
(Modular Earth Submodel System). MADE3 builds on the 
predecessor aerosol submodels MADE and MADE-in. Its 
main new features are the explicit representation of coarse 
mode particle interactions both with other particles and with 
condensable gases, and the inclusion of hydrochloric acid 
(HC1) / chloride (Cl) partitioning between the gas and con- 
densed phases. The aerosol size distribution is represented 
in the new submodel as a superposition of nine lognormal 
modes: one for fully soluble particles, one for insoluble par- 
ticles, and one for mixed particles in each of three size ranges 
(Aitken, accumulation, and coarse mode size ranges). 

In order to assess the performance of MADE3 we compare 
it to its predecessor MADE and to the much more detailed 
particle-resolved aerosol model PartMC-MOSAIC in a box 
model simulation of an idealised marine boundary layer test 
case. MADE3 and MADE results are very similar, except in 
the coarse mode, where the aerosol is dominated by sea spray 
particles. Cl is reduced in MADE3 with respect to MADE 
due to the HC1 / Cl partitioning that leads to Cl removal 
from the sea spray aerosol in our test case. Additionally, 
the aerosol nitrate concentration is higher in MADE3 due 
to the condensation of nitric acid on coarse mode particles. 
MADE3 and PartMC-MOSAIC show substantial differences 


in the fine particle size distributions (sizes < 2 pm) that could 
be relevant when simulating climate effects on a global 
scale. Nevertheless, the agreement between MADE3 and 
PartMC-MOSAIC is very good when it comes to coarse 
particle size distributions (sizes > 2 pm), and also in terms 
of aerosol composition. Considering these results and the 
well-established ability of MADE in reproducing observed 
aerosol loadings and composition, MADE3 seems suitable 
for application within a global model. 


1 Introduction 

Aerosol particles affect the energy balance of the Earth both 

directly by scattering or absorbing radiation and indirectly 
by acting as cloud condensation nuclei. The present-day 

net effect of these processes is probably a negative radia- 
tive forcing (RF) with respect to preindustrial times (e.g. 
Forster et al., 2007; Bellouin et al., 2013; Naik et ah, 2013). 
Hence, the concurrent positive forcing by long-lived green- 
house gases may be partly offset by the aerosol forcing. 

Global model simulations have indicated that the emis- 

sions from ocean ship traffic may be one of the largest con- 
tributors to the anthropogenic aerosol forcing (Lauer et ah, 
2007, 2009; Righi et ah, 2011, 2013; Olivie et ah, 2012; 
Peters et ah, 2012, 2013). That contribution is mainly caused 
by the sulfur in the ship exhaust plumes that leads to the 
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formation of aerosol sulfate (SO 4 ). Both nucleation of new 
particles and condensation of sulfuric acid (H 2 SO 4 ) vapour 
on emitted as well as on background particles contribute to 
SO 4 formation. 

The International Maritime Organization (IMO) has set 
limits on allowed shipping fuel sulfur content that will be fur- 
ther tightened in the future (IMO, 2011) in order to improve 
air quality in port cities and along coasts. On the one hand, 
the sulfur reduction leads to a decrease of aerosol SO 4 con- 
centrations (Lauer et al., 2009; Righi et ah, 2011; Johansson 
et ah, 2013). On the other hand, the reduced SO 4 concen- 
trations allow more aerosol nitrate (NO 3 ) formation by con- 
densation of nitric acid (HNO 3 ). Increased NO 3 content was 
shown to make up for a substantial fraction of the SO 4 re- 
duction (Lauer et ah, 2009; Bellouin et ah, 2011; Righi et ah, 
2011). The aerosol RF and total particulate mass concentra- 
tions may therefore not be reduced as much as the SO 4 con- 
centrations. 

A number of measurements suggest that aerosol NO 3 may 
primarily partition to the coarse mode both under clean ma- 
rine conditions and when marine aerosol is affected by an- 
thropogenic pollution (Kerminen et ah, 1997; Hara et ah, 
1999; Yeatman et ah, 2001; Cavalli et ah, 2004; Nolte et ah, 
2008). However, in the assessment of the ship emissions’ ef- 
fects on climate by Peters et ah (2012, 2013), NO 3 forma- 
tion was not included at all, and the low-sulfur shipping fuel 
studies by Lauer et ah (2009) and Righi et ah (2011) did not 
include interactions of condensable gases with coarse mode 
particles (except for water vapour). These deficiencies may 
have led to errors in the quantification of low-sulfur fuel ef- 
fects. 

To improve on the previous investigations, we devel- 
oped the new aerosol submodel MADE3 for use within 
the MESSy framework (Jockel et ah, 2005, 2010). Via 
this framework, MADE3 can be coupled to an atmospheric 
chemistry scheme and to the general circulation model 
ECHAM5 (Roeckner et ah, 2006), which together form the 
atmospheric chemistry general circulation model (AC-GCM) 
EMAC (ECHAM5/MESSy2 atmospheric chemistry model). 
MADE3 is based on MADE (Ackermann et ah, 1998; Lauer 
et ah, 2005) and MADE-in (Aquila et ah, 2011), but for the 
first time includes interactions of coarse mode particles with 
both condensable gases and other particles. To enable the fur- 
ther extension we restructured and improved major parts of 
the submodel code. Our main motivation for the develop- 
ment of MADE3 is a reassessment of the aerosol perturba- 
tions caused by global ship traffic under current and future 
scenarios for fuel sulfur content. 

Some essential processes, e.g. particle transport and depo- 
sition, are not included in MADE3 because they are treated 
by other submodels within the MESSy framework. It is 
therefore not feasible to evaluate a stand-alone (box model) 
setup of MADE3 by comparison with measured data. In- 
stead, we test here the algorithms used in MADE3 for solv- 
ing the aerosol microphysics equations. In order to assess 
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improvements, strengths, and weaknesses of the new sub- 
model, we compare it to its predecessor MADE and to 
the particle-resolved stand-alone aerosol model PartMC- 
MOSAIC (Riemer et ah, 2009; Zaveri et ah, 2008) in a box 
model application. For that purpose we define a marine back- 
ground setup with added emissions representative of heavy 
ship traffic. We use MADE for comparison because pre- 
vious studies on the shipping effect were carried out with 
this submodel (Lauer et ah, 2007, 2009; Righi et ah, 2011, 
2013). PartMC-MOSAIC is regarded as a reference to test 
how well MADE3 performs as a solver for the aerosol dy- 
namics equation. Implementation and evaluation of MADE3 
within EMAC will be the subject of a follow-up study. 

This paper is organised as follows: in Sect. 2 we first de- 
scribe MADE3 in detail, then briefly state the most important 
improvements with respect to its predecessor MADE, and fi- 
nally summarise the main features of PartMC-MOSAIC with 
a focus on the differences to MADE3. Section 3 contains the 
definition of our marine boundary layer test case. We report 
and discuss the results of the simulations with this setup in 
Sect. 4. Finally, in Sect. 5, we summarise our findings and 
present the conclusions. 

2 Model description 

We use the term “aerosol (sub)modef’ here to refer to the 
computer code used to solve the aerosol dynamics equation. 
Generally speaking, aerosol dynamics includes emissions, 
gas-to-particle conversion, transport, physical and chemical 
processing, and deposition of particles. Throughout this pa- 
per we focus on internal processes, i.e. those that are ac- 
tually calculated by the aerosol submodel MADE3. These 
processes include gas-particle partitioning of semi-volatile 
species, condensation of nonvolatile H 2 SO 4 , formation of 
secondary organics, new particle formation by nucleation, 
and particle coagulation. We add emissions here in order to 
test the reaction of the internal processes to external perturba- 
tions. Transport and loss processes are excluded because they 
are treated by other submodels within the MESSy frame- 
work. 

The aerosol dynamics equation thus takes on the following 
general form (nomenclature based on Riemer et al., 2009): 
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Each aerosol particle is described by a vector n composed 
of the masses fi a of the species a = 1, . . . , A. The number 
distribution of these particles (defined in the A-dimensional 
“species space”) is given by 


t) 


d A N*(fi,t ) 
9/u,i 9/X2 - - - 3 /xa ’ 


( 2 ) 


where N*(fi,t) is the cumulative number concentration of 
particles containing less than the mass /r, fl of species a. 
The two coagulation terms (i.e. gain and loss) are calcu- 
lated based on the collision probability of par- 

ticles described by the vectors fi i and /t 2 - The per-particle 
flux of condensable gases between the gas and the condensed 
phase is /, (fi, g. t ) where i — w stands for water vapour and 
the vector g describes the gas composition. The compo- 
nents of the gas phase vector are the concentrations g, of the 
i = 1, . . . , G different gas species. Condensable gases, i.e. the 

ones considered here, are assumed to be the first i = 1 C 

of these species and correspond to the first a = 1 C 

aerosol species. The species C + 1 is assumed to be water 
vapour or liquid water, respectively. The number distribution 
production rates n* uc (fi,t) and describe the ad- 

dition of new particles by nucleation and emission, respec- 
tively. 


2.1 MADE3 


MADE3 is based on MADE-in (Aquila et al., 2011), an 
extension of MADE as described by Lauer et al. (2005). 
The first generation of MADE was developed for applica- 
tion in a regional model (Ackermann et al., 1998). It was 
derived from work by Whitby et al. (1991) and Binkowski 
and Shankar (1995). Subsequently, MADE was adapted for 
global applications and implemented into the general circula- 
tion model ECHAM4 by Lauer et al. (2005), and later trans- 
formed into a submodel (Lauer et al., 2007) for the MESSy 
framework. The second generation submodel MADE-in was 


developed by Aquila et al. (2011) as an extension to the 
MADE version used by Lauer et al. (2007). It was created 
to enable simulation of number concentrations and mixing 
states of particles containing the insoluble components black 
carbon and mineral dust. For the first version of the third gen- 
eration submodel MADE3, we extended the microphysical 
calculations to also take into account coarse mode particles, 
which were formerly regarded as passive. For version 2.0b, 
we also extended the gas-particle partitioning scheme. The 
“b” stands for beta, which we include in the version number 
because MADE3 has not yet been tested as part of the 3-D 
model. Hence, some minor changes may still be required to 
reach a fully operational version 2.0. 

Despite substantial restructuring of the code during the de- 
velopment of MADE3 and the addition of new features, the 
third generation submodel still shares with MADE most of 
its features and the concepts underlying the computer code. 
The following sections are meant to serve as a reference for 
the mathematics on which the new submodel is based. 


2.1.1 Aerosol properties 


Different numerical representations of the number distribu- 
tion function «*(//,, t) are used in aerosol (sub)models, de- 
pending on the available computational resources, i.e. on the 
target application. A very accurate representation is to track 
particles individually, as is done in PartMC-MOSAlC (see 
Sect. 2.3). However, if N p particles are tracked individually, 
the condensation/evaporation terms in Eq. (1), for instance, 
have to be calculated N p times. To adequately represent the 
whole size range of atmospheric particles, N p has to be of 
the order of 10 5 . The computational cost of such an ap- 
proach is prohibitive for application in an AC-GCM such as 
EMAC. Therefore, the number distribution t) is repre- 
sented in MADE3 in a simplified manner, namely using the 
so-called modal approach. Instead of individual particle com- 
ponent masses, the characteristic variable is taken to be par- 
ticle diameter, and the number distribution is approximated 
by a superposition of nine modes, i.e. nine lognormal func- 
tions (In D, t ) in diameter space (k = 1 , . . . , 9 is the mode 
index): 


;z(lnZ), t ) = 


dN(\nD,t) 

9 In D 


= ^2>n(\nD,t) 


k= 1 
9 


[InD— In £>g 


= y\ Nk{t) e 
\/27r lno-jt 


( 3 ) 


where N(\nD, t ) is the cumulative number concentration of 
particles with diameters smaller than D. The tilde (as in D) is 
used to indicate that the diameter was made dimensionless by 
division by 1 pm. Each mode k is described by three param- 
eters, namely the number concentration Nk(t) of particles in 
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Table 1 . Parameters used in MADE3: mode widths, bulk aerosol 
component densities, accommodation coefficients for gas adsorp- 
tion on aerosol particles, and diffusivities of gases in air. Abbrevi- 
ations are as follows: POM for particulate organic matter, BC for 
black carbon, DU for mineral dust, and SOA for secondary organic 
aerosol. 



Symbol 

Value 

Unit 

Mode widths 

Soluble Aitken 

“1 

1.7 


Mixed Aitken 

“2 

1.7 


Insoluble Aitken 

03 

1.7 


Soluble accumulation 

cr 4 

2.0 


Mixed accumulation 

“5 

2.0 


Insoluble accumulation 

“6 

2.0 


Soluble coarse 

°7 

2.2 


Mixed coarse 

“8 

2.2 


Insoluble coarse 

“9 

2.2 


Component densities 

SO 4 

PSO4 

1.8 x 10 3 

kg m -3 

nh 4 

PNH 4 

1.8 x 10 3 

kg m “ 3 

N0 3 

PNO3 

1.8 x 10 3 

kg m -3 

Na 

PNa 

2.2 x 10 3 

kg m -3 

Cl 

PCI 

2.2 x 10 3 

kg m -3 

POM 

PPOM 

1.0 x 10 3 

kg m -3 

BC 

pbc 

2.2 x 10 3 

kg m -3 

DU 

PDU 

2.5 x 10 3 

kg m -3 

h 2 o 

ph 2 o 

1.0 x 10 3 

kg m -3 

Accommodation coefficients 

H 2 so 4 

“h 2 so 4 

1.0 


nh 3 

“nh 3 

0.1 


hno 3 

“hno 3 

0.1 


HC1 

“HC 1 

0.1 


SOA 

“SOA 

1.0 


Gas diffusivities 

H 2 so 4 

a h 2 so 4 

0.09 

2 —1 
crrr s 

nh 3 

a NH 3 

0.1 

2 — 1 
cm s 

hno 3 

a hno 3 

0.1 

cm s 

HC1 

a HC 1 

0.1 

cm s 

SOA 

a soa 

0.05 

cm s 


that mode, its median particle diameter (geometric mean di- 
ameter) Dg k(t), and its width (geometric standard deviation) 
< 7k . Note that we do not write out time dependencies explic- 
itly in the remainder of this paper. Mode widths are fixed 
in MADE3 (as was the case for its predecessors) in order 
to further reduce the computational burden of the submodel. 
The values of o> are listed in Table 1. All particles in one 
mode k are assumed to have the same composition, i.e. the 
mass fractions c a ^/Yl^-\ c s,k are the same for these parti- 
cles. The symbol c a zk (c s , k ) denotes the mass concentration 
of species a (j) in mode k per unit volume of air. 


Information about particle composition is lost in the size 
distribution representation as given by Eq. (3). Hence, we 
also track the mass concentrations c a ± as described in the 
next subsection. In MADE3 the number of equations to be 
solved is thus 9 x (1 + A), i.e. one equation for the number 
concentration ( A'/.) per mode and A equations for the differ- 
ent aerosol component species ( c a ,k ) per mode. With A — 9 
species only 90 equations are required in MADE3 instead of 
N p ~ 10 5 in PartMC-MOSAIC to solve the aerosol dynam- 
ics equation (Eq. 1). The median diameter D gk of mode k 
can be derived from the component mass concentrations c ak 
in that mode under the assumption of spherical particles ( p a 
is the density of species a, see Table 1): 


D, 


gA 


_( 6V k 

\*N k 


( 4 ) 


with 


Vk = J2 — ( 5 ) 

being the particle volume concentration of mode k. 

Integral moments of the lognormal distribution are often 
used in the internal MADE3 computations to facilitate the 
calculations described in the following subsection. The / th 
moment of mode k is defined as 


M j k = J D' ■ 1 pm 7 • n k (\n D)dln D 

— OO 

= N k .(D g , k y e %^ 2 . 


( 6 ) 


The 0th moment is a mode’s number concentration ( N k ), 
the second moment is related to its particle surface area con- 
centration, and the third moment is related to the mode’s par- 
ticle volume concentration ( V k = tt ■ Mi,^ k /6, see also Eq. 5). 
These moments are used, for instance, to calculate coagula- 
tion and condensation rates (Whitby et al., 1991; Lauer et al., 
2005; Aquila et al., 2011). 

The nine modes in MADE3 (Fig. 1) are representative of 
three size ranges: the Aitken (tens of nanometres), accumula- 
tion (hundreds of nanometres), and coarse modes (microme- 
tres). In each size range MADE3 includes one mode of fully 
soluble particles, one mode of insoluble particles and one 
mixed mode. The choice of this set of aerosol modes follows 
the ideas presented by Aquila et al. (2011), now extended to 
cover also the coarse mode size range. 

Particles in MADE3 consist of up to nine different com- 
ponents, sometimes also called tracers, as they may repre- 
sent more than one chemical species (Fig. 1, Table 1): sulfate 
(SO 4 ), ammonium (NH 4 ), nitrate (NO 3 ), a tracer that con- 
tains sea spray components other than chloride (named Na), 
chloride (Cl), particulate organic matter (POM), black car- 
bon (BC), mineral dust (DU), and water (H 2 0). Note that we 
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free particles 



Internally mixed 
particles with 
a BC/DU core 



Externally 
mixed BC and 
DU particles 



Aitken mode mode coarse mode 


Figure 1 . Schematic illustration of the MADE3 modes and aerosol composition. Each size range (Aitken, accumulation, and coarse mode 
size range) is represented by three modes: soluble (dotted lines), insoluble (dashed lines), and mixed particles (solid lines). Red lines indicate 
the new modes in MADE3 with respect to its predecessor MADE-in. Small “pie charts” show possible Aitken mode particle compositions, 
larger “pie charts” stand for accumulation and coarse mode particles. Note that the MADE3 code in principle allows all components in each 
of the modes, but some components do not appear in significant amounts in certain types of particles in the real atmosphere (e.g. mineral 
dust in Aitken mode particles). Abbreviations are as follows: POM for particulate organic matter, BC for black carbon, and DU for mineral 
dust. 


do not add charges to the tracer names here because we do 
not distinguish between different oxidation and phase states 
of the aerosol components. For instance, the SO 4 tracer in- 
cludes SO^ - , HSO^ , and liquid H 2 SO 4 , as well as the sulfate 
in (NH 4 ) 2 SC >4 and other crystalline salts. 

2.1.2 Aerosol processes 

Particle composition, number concentration, and size distri- 
bution undergo changes during the atmospheric processing 
of the aerosol. MADE3 calculates the evolution of the parti- 
cle population due to three main processes: ( 1 ) gas-particle 
partitioning of semi-volatile species and water, which adds 
or removes particulate mass depending on environmental pa- 
rameters (such as temperature or relative humidity); ( 2 ) gas- 
to-particle conversion of low-volatility species by conden- 
sation on pre-existing particles or in situ formation of new 
particles; (3) mass transfer between particles by coagulation, 
which concurrently reduces particle number concentration. 
Employing an operator splitting approach, MADE3 first cal- 
culates compositional changes due to gas-particle partition- 
ing alone. Subsequently, condensation/new particle forma- 
tion and coagulation are treated simultaneously. 

Gas-particle partitioning 

Gas-particle partitioning of semi-volatile trace constituents 
(NH 3 /NH 4 , HNO 3 /NO 3 , HC1/C1) and water is calculated 
in MADE3 using the thermodynamic equilibrium model 
EQSAM (EQuilibrium Simplified Aerosol Model; Metzger 


et al., 2002, 2006). A sequential treatment of the three size 
ranges is applied in analogy to the procedure adopted by 
Aquila et al. (2011). The assumption of equilibrium between 
the gas and condensed phases is well justified for fine par- 
ticles because they equilibrate on timescales of the order of 
seconds up to a few minutes (e.g. Meng and Seinfeld, 1996), 
i.e. well within the typical time step used for aerosol studies 
with EMAC (~ 30 min in T42 resolution, i.e. ~ 2.8° x 2.8°). 

For large particles (in the size range of micrometres and 
greater), however, gas diffusion may be too slow to en- 
able equilibration within this time frame (e.g. Wexler and 
Seinfeld, 1990, 1992; Meng and Seinfeld, 1996). Flence, 
the equilibrium assumption may introduce substantial er- 
rors (Moya et al., 2002; Koo et al., 2003; Feng and Penner, 
2007; Athanasopoulou et al., 2008), but fully dynamical cal- 
culations of the involved fluxes are infeasible in long-term 
simulations with AC-GCMs. Several solutions to this prob- 
lem have been proposed (Capaldo et al., 2000; Pilinis et al., 
2000; Jacobson, 2005; Zaveri et al., 2008). However, most 
of them are still computationally too expensive for appli- 
cation in MADE3 within EMAC, and the solution by Za- 
veri et al. (2008) would require a complete revision of the 
chemistry scheme that was used in previous simulations with 
MADE (e.g. Lauer et al., 2007; Aquila et al., 2011; Righi 
et al., 2011, 2013). Consequently, we account for possible 
non-equilibrium effects by limiting the gas-particle fluxes 
involving coarse mode particles in a similar manner as de- 
scribed by Pringle et al. (2010a); Pringle et al. (2010b). We 
calculate the maximum possible diffusion fluxes of the semi- 
volatile gases (using equations equivalent to Eqs. ( 8 ) and (9) 
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with gx, s = 0 , see below) and perform an equilibrium calcu- 
lation with EQSAM. If the fluxes required to reach this equi- 
librium surpass the maximum fluxes, we limit the amount of 
material that condenses during one time step accordingly. 


Condensation of H 2 SO 4 and organic vapours 


Due to its very low equilibrium vapour pressure, we assume 
that all H 2 SO 4 is transferred from the gas phase to the aerosol 
phase during each time step. Depending on the magnitude of 
the condensation flux, this transfer can occur via condensa- 
tion alone, or via condensation and new particle formation 
(see below). To determine the amount of H 2 SO 4 that can 
condense during one time step, we calculate the condensation 
flux explicitly. This is also necessary to obtain the proper dis- 
tribution of the condensate among the differently sized parti- 
cles. An equilibrium assumption does not yield this distribu- 
tion. While potential errors in the distribution of condensing 
material among differently sized particles will be corrected 
by re-evaporation in the case of the semi-volatile species, 
this is not possible for H 2 SO 4 , since our assumption that all 
H 2 SO 4 is transferred to the aerosol phase means that it can- 
not re-evaporate in our model. The total condensation flux of 
H 2 SO 4 is the sum of the rates of change of mass concentra- 
tions CH 2 so 4 ,/t f°r all modes (k = 1, . . . , 9): 


/ dcH 2 S0 4 T 

V dr 


cond 


r dm p (D) 

J dr 


0 


n k (D)dD, 


(7) 


where dm p (D)/dr is the rate of change of mass for 
an individual particle of diameter D and n k (D)dD = 
n k (\nb)d\nD. Note that with dlnD/dZ) = 1/Z) one ob- 
tains n k ( In D)/D as the functional form for n k (D), where 
n k (lnD) has to be inserted as given in Eq. (3). The rate of 
particle mass change depends on the ratio of the particle di- 
ameter D and the mean free path X of H 2 SO 4 molecules in 
the gas phase. As atmospheric aerosol particles span a large 
range of sizes, two limiting cases have to be considered. In 
the continuum regime, i.e. for D ~S> X, one obtains (see Sein- 
feld and Pandis, 2006) 


{ dnip(D) 
V dr 


cont 

= 27tDAh 2 so 4 C?h 2 so 4 ,00 


gH 2 S0 4 ,s) 


( 8 ) 


with the gas phase difflisivity Ah 2 so 4 (see Table 1). The ad- 
ditional indices “ 00 ” and “s” specify gas phase concentra- 
tions far away from the particle surface and directly above it, 
respectively. The corresponding expression for the kinetic, or 
free molecular, regime where D <^X is (Seinfeld and Pandis, 
2006) 

^ d/n p (D) y~ ree = 

7 r D 2 

— «H 2 S0 4 ®H 2 S0 4 (<?H 2 S0 4 ,oo - gH 2 S0 4 ,s)- (9) 


Here, an,so 4 is the accommodation coefficient (see Ta- 
ble 1) and ®h 2 so 4 the thermal speed of the H 2 SO 4 molecules: 


ft) H 2 S0 4 = 


I 8 RT 
xMh 2 sc>4 


( 10 ) 


where R is the universal gas constant, T the absolute tem- 
perature, and Mh 2 so 4 the molar mass of H 2 SO 4 . Note that 
a wide range of values for aH,so 4 h as been derived from 
measurements (0.02-0.79; e.g. Van Dingenen and Raes, 
1991; Kerminen and Wexler, 1995; Jefferson et al., 1997; 
Bardouki et al., 2003) and used in models (0.1-1; e.g. 
Capaldo et al., 2000; Vignati et al., 2004; Zaveri et al., 
2008; Mann et al., 2010; Kajino et al., 2012). Here, we use 
«h 2 S0 4 = 1 as in former generations of the submodel, which 
was also found to be in agreement with field measurements 
by Eisele and Tanner (1993). The value is the same for all 
modes, i.e. condensation is treated in the same way, regard- 
less of whether particles contain insoluble material or not. 
We set #h 2 so 4 ,s = 0 due to the very low equilibrium vapour 
pressure of H 2 SO 4 . For gH 2 so 4 ,oo we use the solution to 
the ordinary differential equation that describes the tempo- 
ral evolution of the gas phase H 2 SO 4 concentration gH 2 S0 4 : 


dgH 2 so 4 

dr 


-Ph 2 S0 4 ~ £h 2 S0 4 • £h 2 so 4 - 


( 11 ) 


Here, / j h 2 so 4 is the production rate of gaseous H 2 SO 4 and 
f n 2 so 4 is the sum of the integrals as given in Eq. (7) for all 
nine modes with the factor £h 2 S0 4 ,oo removed, i.e. the overall 
loss coefficient due to condensation: 


Th 2 so 4 = 


v-9 /dc H ,so 4 .r\ 

\ / cond 


gH 2 S0 4 ,oo 


( 12 ) 


Production of gaseous H 2 SO 4 and condensation on the 
particles are thus considered in parallel. 

The integral in Eq. (7) can be evaluated analytically with 
both expressions for dm p (D)/dr (Whitby et al., 1991). How- 
ever, there is also a transition regime, where D and X are of 
the same order of magnitude. We therefore apply the method 
of Binkowski and Shankar (1995), using half the harmonic 
mean of the two integrals as the condensation rate of H 2 SO 4 
on mode k. For more details on the condensation calculations 
see the description in Appendix A of Aquila et al. (2011) or 
the original work by Whitby et al. (1991). 

Secondary organic aerosol (SOA) formation from con- 
densing organic vapours is treated in the same simplified 
manner in MADE3 as in MADE-in (Aquila et al., 2011). 
Using an externally supplied mass formation rate of SOA 
( P SO a), we apply a similar procedure as outlined above for 
H 2 SO 4 condensation. The SOA production rate is multiplied 
by the time step duration A t to obtain g soa.oo for analo- 
gous expressions to Eqs. ( 8 ) and (9). The near-surface gas 
phase concentration is again set to zero (gsoA.s = 0, neglect- 
ing semi-volatile organic species because organic gas phase 
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chemistry is not considered) so that the integral in Eq. (7) can 
be evaluated just as for H 2 SO 4 . 

New particle formation 

Nucleation of new particles from H 2 SO 4 and H 2 O is calcu- 
lated in MADE3 after solving the production-condensation 
equation (Eq. 11). This approach corresponds to method 2C 
as discussed by Wan et al. (2013). In terms of the nucleation 
sink for gaseous H 2 SO 4 , they showed this method to be the 
best-performing among sequential methods for solving the 
full H 2 SO 4 equation, i.e. Eq. (11) plus a nucleation loss term. 
To calculate the nucleation rate, (dfVi/dr) nuc , we apply the 
parameterisation by Vehkamaki et al. (2002, 2013) that is 
based on temperature, relative humidity (RH), and H 2 SO 4 
concentration. Following Binkowski and Roselle (2003), we 
account for rapid growth of the freshly nucleated particles 
to detectable sizes by assuming a monodisperse size distri- 
bution with D — 3.5 mn upon formation. The H 2 SO 4 frac- 
tion of these particles is calculated from the ambient RH as 
described in Binkowski and Roselle (2003), based on mea- 
surements by Nair and Vohra (1975). Subsequently, number 
and mass concentrations of the nucleated particles are added 
to the soluble Aitken mode ( k = 1). New particle formation 
from organic precursor gases is not considered in MADE3. 

Coagulation 


The parameters pi and p m in Eq. (14) stand for the densi- 
ties of particles in modes / and m, respectively. The coeffi- 
cients af m are calculated as follows: 

as) 

Here, as in Eq. (14), the Kronecker symbol S x y has been 
used. Its value is 8 Xi y =1 if x = y, and S Xt y = 0 otherwise. 
The matrix elements r /,„ are used for the assignment of num- 
ber and mass concentrations of coagulated particles to a tar- 
get mode, depending on the modes of origin, / and in. The 
rules for this assignment, i.e. the values of the matrix ele- 
ments Tim, follow Aquila et al. (2011): 

- intramodal coagulation produces particles in the same 
mode (t// = /); 

- intermodal coagulation produces particles in the size 
range of the larger mode; 

- the exact target mode for intermodal coagulation de- 
pends on the mass fraction x of soluble material and 
water in the final particles: 

x = 1 : soluble mode, 

0.1 < x < 1: mixed mode, 
x < 0.1: insoluble mode. 


Similar to the condensation treatment, coagulation calcula- 
tions in MADE3 are also performed by mode. Number and 
mass changes are calculated separately: 



(13) 


CXJ cx, 

// 

0 0 


P(D\, D 2 )ni(D\)n m (D 2 )dDi dD 2 


ic a .k 
~“d T 


9 9 


= r££ L 

coag u /= 1 m =l L 


{h,T lm ~ S k ,l) 


Ca,l 


E A 

j= 1 c -s,/ 


OO OO 

■PI- j J(D l ) 3 p(D l ,D 2 )ni(Di)n m (D 2 )dDidD2 
0 0 


+ ( &k,Ti m &k,m) ' 


E a 

s= 1 C - 5 , 


(14) 


■Pm- j J(D 2 ) 3 P(D l ,D 2 )n l (D l )n m (D 2 )dD l dD2 
00 


For instance, particles that result from intermodal coagula- 
tion of particles from modes l = 1 (soluble Aitken mode, or 
“ks” in Table 2) and m — 4 (soluble accumulation mode, or 
“as”) are assigned to mode 744 = 4 (“as”). Hence, af 4 = 0, 
which means that this process does not add to the particle 
number in mode k — 4 (“as”). It does, however, add mass 
from mode / = 1 (“ks”) to mode k — 4 (“as”). This is re- 
flected in the parentheses with the Kronecker symbols in 
Eq. (14): the first pair of parentheses evaluates to one, the 
second pair to zero. In case of intramodal coagulation, i.e. if 
l = k and m = k , the value of the coefficient in Eq. (13) is 
a kk = —0.5. It is negative because one particle per such event 
is lost, but the factor is only —0.5 because of the double inte- 
gration over the same mode. For intramodal coagulation, all 
the Kronecker symbols in Eq. (14) evaluate to one, so that 
all summands are zero, and no mass is added to, or removed 
from mode k. 

In total, Eq. (13) and Eq. (14) include 45 and 90 sum- 
mands for each mode k, respectively, but many of them are 
zeros. For example, coagulation losses and gains in the sol- 
uble coarse mode (“cs” in Table 2) are described by seven 
non-zero terms in the number equation (Eq. 13), and eight 
non-zero terms in the mass equation (Eq. 14). 

For the Brownian coagulation kernel P(D\,D 2 ) we use 
the approximate formulations developed by Whitby et al. 
(1991) that can be integrated analytically. Two different ex- 
pressions are required again, depending on the size of the 
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Table 2. Matrix elements r/ m (table cells) for assignment of particles that result from coagulation of particles from mode Z (row) with particles 
from mode m (column). The mode naming convention is “k”, “a”, “c” as a first letter to specify the Aitken, accumulation, or coarse mode, 
respectively, and “s”, “m”, “i” as a second letter for the soluble, mixed, or insoluble mode, respectively. Thus, the soluble coarse mode, for 
instance, is named “cs”. Corresponding numbers are as follows: ks = 1, km = 2, ki = 3, as = 4, am = 5, ai = 6, cs = 7, cm = 8, and ci = 9. 



ks 

km 

ki 

as 

am 

ai 

cs 

cm 

ci 

ks 

ks 

km 

km/ki 

as 

am 

am/ai 

cs 

cm 

ci 

km 

km 

km 

km/ki 

am 

am 

am/ai 

cm 

cm 

ci 

ki 

km/ki 

km/ki 

ki 

am/ki 

am/ki 

ai 

cm 

cm 

ci 

as 

as 

am 

am/ki 

as 

am 

am/ai 

cs 

cm 

cm/ci 

am 

am 

am 

am/ki 

am 

am 

am/ai 

cm 

cm 

cm/ci 

ai 

am/ai 

am/ai 

ai 

am/ai 

am/ai 

ai 

cm/ai 

cm/ai 

ci 

cs 

CS 

cm 

cm 

CS 

cm 

cm/ai 

CS 

cm 

cm/ci 

cm 

cm 

cm 

cm 

cm 

cm 

cm/ai 

cm 

cm 

cm/ci 

ci 

ci 

ci 

ci 

cm/ci 

cm/ci 

ci 

cm/ci 

cm/ci 

ci 


particles. For the continuum regime the function is given as 
(Whitby et al., 1991) 

p con \D l ,D 2 ) = 


2 k B T 

2 + 2 A. A 
(' L 

( 1 

1 02 

3v 

+ 21A 

Ui 
>1 ^ 

CDi) 2 
+ ° 2 + 

K d 2 1 (d 2 ) 2 ) 

+ D, + 


D 2 . 


(16) 


with Boltzmann’s constant k\i, atmospheric dynamic viscos- 
ity v, and the constant A = 1.246 that accounts for the re- 
duced drag on small particles. Atmospheric viscosity is cal- 
culated from temperature (T): 


v = B- 


3 

T2 

T + S 


(17) 


with B= 1.458 x 10 -6 kg m -1 s” 1 K -0 5 and S= 110.4 K. 
The mean free path depends on both temperature and pres- 
sure: 


k = A- 


PqT 
Top ' 


(18) 


where A = 6.6328 x 10 -8 m, po = 101 325 Pa, and 7o = 
288. 15 K. For the free molecular regime, the coagulation ker- 
nel becomes (Whitby et al., 1991) 


f} bee (D\, D 2 ) = 


6k B T 

P 1 +P 2 



(Pi ) 2 

(Di) 3 / 2 


(19) 


QPi ) 2 

(Dt) 3 ! 2 


)_R}_ 

Sd; 



where p\ (p 2 ) is the density of the particle with diameter I)\ 
(D 2 ). Note that a correction factor is required for the inte- 
grals in Eqs. (13) and (14) if this kernel approximation for 
the free molecular regime is used. It is set constant at 0.8 for 
unimodal and 0.9 for bimodal coagulation in MADE3. Af- 
ter evaluation of the coagulation integrals for both regimes 


(i.e. with /i com (A)| , D 2 ) and /3 frec (/9| , D 2 ), respectively) the 
halved harmonic means of the resulting rates are used to 
redistribute mass and numbers among the modes (see Ap- 
pendix B in Aquila et al., 2011). 

Renaming 

As particles grow by condensation and coagulation, the 
Aitken modes may grow into the size range of the accumula- 
tion modes. In order to avoid mode merging, i.e. to keep the 
modes approximately within their assigned size ranges, we 
apply a procedure that is termed renaming (Binkowski and 
Roselle, 2003). One of two criteria must be met in MADE3 
to trigger renaming within a time step. Either the volume 
growth rate of the Aitken mode must be larger than that of 
the corresponding accumulation mode, or the median diam- 
eter of the Aitken mode must exceed 30 mn and its number 
concentration must be greater than that of the corresponding 
accumulation mode. In such a case the number concentration 
of particles greater than the intersection diameter of the two 
number size distributions is shifted from the Aitken to the 
corresponding accumulation mode. The associated mass con- 
centration is also transferred. Renaming is performed only 
between modes of the same particle type, i.e. either between 
the two soluble modes, or between the two insoluble modes, 
or between the two mixed modes. Note that we do not rename 
particles from the accumulation to the coarse modes because 
their diameters are changed much less by condensation and 
coagulation than those of the Aitken mode particles. 

Aging of insoluble particles 

The aerosol processing in the atmosphere is also termed ag- 
ing. For insoluble particles this term often refers to the ac- 
quisition of a coating of soluble components that transforms 
them from an initially hydrophobic state to a hydrophilic one. 
In MADE3 this transformation is realised by transfer of num- 
ber and mass concentrations from the insoluble modes to 
the mixed modes in analogy to the procedure described by 
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Aquila et al. (2011). Following Aquila et al. (2011), we use 
a threshold mass concentration fraction of 10% of soluble 
material and water in an insoluble mode to trigger this trans- 
fer. 

2.2 MADE3 vs. MADE 

MADE is used here in the version described by Lauer 
et al. (2005). Although the code underwent major restruc- 
turing and was expanded for the development of MADE3 via 
MADE-in (Aquila et al., 2011), the new submodel still shares 
with MADE the computational approaches to aerosol size 
distribution representation, gas-particle partitioning, H2SO4 
and SOA condensation, new particle formation, and coagu- 
lation. Here we therefore point out only the major physically 
motivated changes in MADE3 with respect to MADE. 

The most obvious difference is in the number of modes, 
which increased from three (MADE) to nine (MADE3). The 
representation of aerosol particles by three modes per size 
range allows us to model both internally mixed particles 
and externally mixed particle populations (see Aquila et al., 
2011 ). 

Furthermore, while particles in the MADE coarse mode 
are considered passive (only water uptake by coarse mode 
particles is included), they interact with both other particles 
and condensable trace gases in MADE3. Coarse mode par- 
ticle composition as well as effects of the coarse mode on 
fine particles and the gas phase can therefore be resolved in 
much more detail with MADE3 than what is feasible with 
MADE. Note also that coarse mode particles in MADE are 
composed of sea spray, mineral dust, and water only, while 
MADE3 allows all aerosol components to be present in the 
coarse modes. 

In combination with the larger number of modes, the 
newly introduced interactions of coarse mode particles also 
entail a larger number of different possible coagulation path- 
ways: 45 in MADE3 vs. 2 in MADE. The calculations to 
determine target modes based on the soluble mass fraction of 
the coagulated particles is not necessary in MADE, while it 
is required for 14 of the coagulation pathways in MADE3. 

In addition. Cl is considered as a separate species in 
MADE3, whereas all sea spray components are lumped into 
one tracer in MADE. The explicit Cl representation en- 
ables the calculation of gas-aerosol partitioning for HC1 by 
EQSAM, which is not considered in MADE, but is required 
for accurate modelling of processes in the marine boundary 
layer and in coastal areas. 

With the larger number of modes, the inclusion of the 
coarse mode particle interactions, and the inclusion of the Cl 
tracer, the number of aerosol species tracers increased from 
18 in MADE to 81 in MADE3. 


2.3 PartMC-MOSAIC 

PartMC-MOSAIC is a stochastic particle-resolved aerosol 
model that consists of the microphysics code PartMC 
(Particle-resolved Monte Carlo model, Riemer et al., 2009) 
and the gas and condensed phase chemical solver MOSAIC 
(MOdel for Simulating Aerosol Interactions and Chemistry, 
Zaveri et al., 2008). The PartMC-MOSAIC version used for 
the present study (v. 2.2.1) corresponds to the detailed de- 
scription in Tian et al. (2014), so that we only give a brief 
summary of the relevant features here. 

The model solves the aerosol dynamics equation (Eq. 1) 
on a per-particle basis. While the size distribution is con- 
strained in MADE3 by the assumption of lognormal modes, 
it can freely evolve in PartMC-MOSAIC, where t ) is 

represented by a finite number /V p of computational particles 
with discrete sizes. For the present study we chose N ? ~ 10 5 , 
as was done in previous applications of PartMC-MOSAIC 
(e.g. Riemer et al., 2009; Tian et al., 2014). In order to cap- 
ture the large range of possible sizes and concentrations, one 
computational particle can represent a larger number of real 
particles (DeVille et al., 2011). The number and mass weight- 
ing of these computational particles for the microphysics and 
chemistry calculations is performed automatically. As parti- 
cles are constantly emitted but loss processes are not consid- 
ered here (except for coagulation), half of the particles are 
randomly picked out and discarded whenever the number of 
computational particles exceeds twice its initial value. 

Aerosol composition can be resolved into more separate 
species in PartMC-MOSAIC than in MADE3. Here, we use 
11 tracers: SO4, NH4, NO3, Na, Cl, organic carbon (OC), 
black carbon (BC), calcium (Ca), carbonate (CO3), other in- 
organic material (OIN), and H2O. 

Particle emissions and coagulation are treated stochasti- 
cally in PartMC. Random samples are added at each time 
step such that the number of emitted particles per unit time 
is Poisson distributed around a prescribed continuous mean 
emission rate. The composition and mean size distribution of 
these particles are also prescribed. For coagulation the max- 
imum number of collision events during the time step is es- 
timated and a corresponding number of candidate particle 
pairs is randomly selected. Subsequently, an accept-reject 
procedure is applied to determine whether these pairs actu- 
ally coagulate. The probability for acceptance is based on the 
Brownian coagulation kernel. 

Condensation of H2SO4 and gas-particle partitioning of 
semi-volatile gases is dynamically calculated by the deter- 
ministic model MOSAIC. This is in contrast to the equilib- 
rium assumption in MADE3, so that no special treatment of 
large particles is required here. Note that, besides the dif- 
ferent approaches to aerosol microphysics, the use of differ- 
ent codes for the thermodynamic calculations (EQSAM in 
MADE3 vs. MOSAIC in PartMC-MOSAIC) can be a major 
driver of differences in simulation results. 
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The validity of the PartMC microphysics routines was 
demonstrated by Riemer et al. (2009) and MOSAIC was 
shown to perform well in comparison to other aerosol chem- 
istry codes (Zaveri et al., 2008), which included more de- 
tails than the treatment in MADE3. The combined model 
(PartMC-MOSAIC) was successfully applied in a recent 
study of a ship plume (Tian et al., 2014). In summary, 
PartMC-MOSAIC is able to capture many more details of 
the aerosol evolution than MADE3 and can therefore serve 
as a reference for our comparison. 

3 Test case setup 

We performed test simulations with MADE3, MADE and 
PartMC-MOSAIC using initial conditions representative of 
the marine background boundary layer. The setup was de- 
signed to mimic the actual target application of MADE3, 
namely its use within the framework of the AC-GCM EMAC. 
As a first application, we plan to use MADE3 for a re- 
assessment of the shipping effect on the global atmospheric 
aerosol, because coarse mode particle interactions with con- 
densable gases and with small particles may play an impor- 
tant role in that context. Therefore, we added black carbon 
(BC) emissions, and prescribed gaseous H 2 SO 4 and HNO 3 
production rates in our test case, thus simulating an episode 
of heavy ship traffic. 

We simulated 24 h of aerosol processing (without trans- 
port and deposition, see Sect. 2) under constant environmen- 
tal conditions, with a constant BC emission rate and constant 
H 2 SO 4 and HNO 3 formation rates. This scenario can be re- 
garded as an idealised representation of a stagnant air mass 
in a shipping corridor. For the time steps we chose typical 
(model-specific) values: 1800 s, i.e. 30 min, in MADE3 and 
MADE, and 1 s in PartMC-MOSAIC. Gas phase chemistry is 
not considered, because we want to focus on the particulate 
phase here. 

Environmental parameters as well as initial gas and 
aerosol concentrations were extracted from a previous multi- 
year EMAC simulation using the MADE3 predecessor 
MADE-in (evaluated in Aquila et al., 2011). For our test case, 
we picked a grid box in the Indian Ocean with T = 286 K, 
p = 1.02 x 10 5 Pa, and RH = 0.771. The initial aerosol state 
(as represented in MADE3) is shown in Fig. 2, and the ini- 
tial number and species mass concentrations per size range, 
i.e. summed over the respective soluble, mixed, and insolu- 
ble modes, are given in Table 3. Initial gas concentrations 
are as follows: [H 2 SO 4 ] = 3.75 x 10 - 7 pg m -3 , [NH 3 ] = 
0.240 pgm- 3 , [HNO 3 ] = 0 , and [HCl] = 0 . 

As MADE-in represents coarse mode particles by only 
one mode, we redistributed the species mass concentrations 
among the MADE3 coarse modes as follows: 


J. C. Kaiser et al.: MADE3 - description and test 

Table 3. Initial number and species mass concentrations in the 
Aitken, accumulation, and coarse mode size ranges (summed over 
the respective soluble, mixed, and insoluble modes). Abbreviations 
are as follows: POM for particulate organic matter, BC for black 
carbon, and DU for mineral dust. 


Aitken mode Accumulation mode Coarse mode 
Number concentrations [m 3 ] 

7.34 x 10 7 3.51 x 10 6 3.44 x 10 6 


Species mass concentrations [|ig nr 


S0 4 

2.37 x 10“ 4 

0.0425 

0 

nh 4 

8.89 x 10“ 5 

2.86 x 10“ 8 

0 

no 3 

0 

0.0161 

0 

Na 

0 

0.161 

6.51 

Cl 

0 

0.200 

8.10 

POM 

4.17 x 10“ 5 

9.31 x 10“ 5 

0 

BC 

1.12 x 10“ 6 

7.24 x 10“ 6 

0 

DU 

0 

2.42 x 10“ 6 

0.00266 

h 2 o 

0 

1.21 

44.1 


- sea spray: 50 % to the soluble mode, 50 % to the mixed 
mode; 

- mineral dust (DU): 50 % to the mixed mode, 50 % to 
the insoluble mode; 

- H 2 O: approximately 50 % each to the soluble mode and 
the mixed mode and 1 .6 x 10 -4 % to the insoluble mode 
(in order to keep the H 2 O mass fraction of the latter 
below the 10 % threshold upon initialisation). 

Other species are not included in the initial coarse mode 
particle composition because MADE-in does not allow sim- 
ulation of other components in the coarse mode. Splitting 
up the MADE-in sea spray tracer, we assigned 45 % of the 
mass concentration to the MADE3 Na tracer and 55 % to the 
MADE3 Cl tracer. This speciation is in accordance with the 
assumptions in EQSAM on sea spray composition. The re- 
distribution of number concentration was derived from the 
mass concentrations in the coarse modes under the assump- 
tion that all three modes should initially have the same me- 
dian diameter. 

Transformation of the initial aerosol state to the MADE 
representation is straightforward: mass and number concen- 
trations from the MADE-in Aitken modes were summed 
up and assigned to the MADE Aitken mode, and the same 
procedure was applied to the accumulation modes. For the 
coarse mode, the MADE-in output could be used without 
modifications. 

In terms of median diameters, number concentrations, 
and mode widths, PartMC-MOSAIC was initialised with the 
same modes as MADE3, translated to a population of indi- 
vidual particles. However, the MADE3 Na and DU tracers 
had to be further speciated for use with PartMC-MOSAIC. 
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Figure 2. Initial aerosol mass size distribution and composition as represented in MADE3. The thick black line represents the mass size 
distribution calculated as the sum of the modes shown in grey (left vertical axis; dotted lines for soluble modes, solid lines for mixed modes, 
dashed lines for insoluble modes). Insoluble Aitken and accumulation mode mass concentrations are initially so small that these modes do not 
appear in the figure, and the curves for the soluble and mixed coarse modes lie on top of each other. The coloured bars show the contributions 
of the individual species (POM = particulate organic matter, BC = black carbon, DU = mineral dust) to the mass concentration (right vertical 
axis) of the respective mode (from left to right: soluble Aitken, mixed Aitken, soluble accumulation, mixed accumulation, soluble coarse, 
mixed coarse, insoluble coarse). Note that while the right vertical axis is logarithmic, the species fractions in the bars add up linearly to the 
total mass concentrations (i.e. the axis only applies to the total mass concentration in each mode, but not to the individual contributions). 
Note further that the three coarse mode bars were artificially spread out along the diameter axis and grey borders corresponding to the line 
styles of the respective modes were added for clarity. 


Following again the sea spray composition assumptions in 
EQSAM, we assigned 69% of the MADE3 Na tracer mass 
concentration to the PartMC-MOSAIC Na tracer, 17% were 
added to the PartMC-MOSAIC SO4 tracer, 3 % to the Ca 
tracer, and 11 % to the other inorganic material (OIN) tracer. 
For the speciation of the MADE3 DU tracer into PartMC- 
MOSAIC tracers we assumed the following mass fractions: 
2 % of Ca, 3 % of CO3, and 95 % of OIN (corresponding to 
5 % CaC03, based on data in Glaccum and Prospero, 1980; 
Kandler et al., 2009; Scheuvens et al., 2013). 

BC emissions are added to the (insoluble) Aitken and ac- 
cumulation modes in MADE (MADE3), and as separate par- 
ticles in PartMC-MOSAIC. Our BC emission flux (see Ta- 
ble 4) is based on the values for ship emissions reported in the 
Lamarque et al. (2010) data set for the year 2000, assuming 
a typical marine boundary layer height of 500 m. We chose 
values from a grid box off the coast of Normandy, France, 
and assumed a bimodal size distribution for the emitted BC 


particles as used in Righi et al. (2013). For the formation rates 
of gaseous H2SO4 and HNO3 (Table 4) we assumed rather 
high values in order to fully explore the effects of the con- 
densation process on the microphysical calculations. These 
rates correspond to a direct conversion of the sulfur diox- 
ide (SO2) and nitrogen oxides (NO x ) emitted by ships (again 
drawn from the Lamarque et al., 2010, data set), i.e. their 
choice implicitly contains the assumption of equilibrium SO2 
and NO x concentrations with respect to ship emissions and 
formation of gaseous H2SO4 and HNO3. 

New particle formation (NPF) is treated differently in 
MADE/MADE3 and PartMC-MOSAIC. Due to the large un- 
certainties associated with the choice of parameterisations, 
consideration of NPF would make the interpretation of the 
simulation results rather difficult. We therefore neglect the 
process here. In a sensitivity experiment we switched on the 
nucleation calculation in MADE and saw a NPF event af- 
ter ~ 3 h of simulated time. We added the number and mass 
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Table 4. Emission and formation rates used in the test case. 


Species 

Mass cone, 
rate of change 
[kg m -3 s -1 ] 

Number cone, 
rate of change 
[m -3 s -1 ] 

Aitken mode BC 

1.9 x 10 “ 16 

2.6 x 10 2 

(Dg = 70 nm, a = 1.45) 
Accumulation mode BC 

5.0 x 10 “ 17 

2.0 

(Do — 260 nm, a = 1.25) 
H 2 SO 4 

1.5 x 10 “ 4 * * * * * * * * * 14 

_ 

HNO 3 

1.7 x 10 “ 14 

- 


size distributions of these nucleated particles as simulated by 
MADE to the initial aerosol state for all three (sub)models. 
Subsequently, we switched NPF off again and re-ran the sim- 
ulations. Differences in the 24 h number and mass size distri- 
butions between these and the original simulations were neg- 
ligible for all (sub)models. Hence, we can assume that the 
atmospheric processing of nucleated particles is adequately 
treated by MADE3, i.e. growth by condensation and removal 
by coagulation with larger particles are properly represented. 

Since some processes are treated by stochastic approaches 
in PartMC-MOSAIC we ran the model ten times and calcu- 
lated the average of the aerosol mass and number concen- 
trations for the ensemble of simulations. This procedure en- 
sures that we do not discuss an “outlier” here and enables us 
to quantify uncertainties. 

4 Results and discussion 

4.1 Comparison of MADE3 and MADE 

4.1.1 Size distributions 

Number and dry mass size distributions in MADE3 and 

MADE at the beginning and at the end of the 24 h simulation 

are plotted vs. dry diameter in Fig. 3. Dry quantities are cal- 
culated from all aerosol components except water. Although 

water constitutes the largest fraction of the aerosol mass (see 

Fig. 4) we chose the dry representation here because the large 

H 2 O mass masks finer features in the size distribution. As 
aerosol water content is diagnosed from the composition of 

the dry aerosol anyway, no essential information is lost from 
the size distributions when neglecting H 2 O here. Deviations 
in the size distributions after 24 h in Fig. 3 are small and can 
be explained by the new features of MADE3 as follows. 

In the fine particle size range we see a bimodal shape of 
the 24 h MADE number size distribution which is not visible 
in the corresponding MADE3 distribution. This difference is 
due to different overlaps of the initial Aitken and accumula- 
tion modes in the two submodels, which leads to a stronger 
convergence of median mode diameters in MADE3 than in 
MADE. For the initial MADE Aitken (accumulation) mode, 
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Size distributions: MADE3 vs. MADE 
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Figure 3. Initial (grey, black) and final (light red, red) size distri- 
butions in MADE3 and MADE vs. dry diameter (i.e. neglecting 
aerosol water). Light colours (grey, light red) are used for num- 
ber distributions (left vertical axis), full colours (black, red) for dry 
mass distributions (right vertical axis). MADE3 output is shown as 
solid lines, MADE output as dashed lines. The inset magnifies the 
coarse mode dry mass size distributions. 


we summed up species mass and number concentrations of 
the initial MADE3 Aitken (accumulation) modes. The as- 
sociated mixing of particles in MADE leads to greater me- 
dian diameters for both modes with respect to the median 
diameters of the soluble Aitken and accumulation modes in 
MADE3. These soluble modes overlap more strongly dur- 
ing the first hours of the simulation than the Aitken and the 
accumulation mode in MADE. Hence, more particles are re- 
named from the Aitken mode to the accumulation mode in 
MADE3 during that time and the accumulation mode is thus 
shifted towards smaller diameters. The renaming stops when 
the number concentration of accumulation mode particles 
surpasses that of the Aitken mode. This happens after about 
14 h and after about 18 h of simulated time in MADE3 and 
MADE, respectively. Subsequently, the Aitken mode parti- 
cles grow towards the accumulation modes by condensation. 
As this growth begins earlier in the MADE3 simulation than 
in the MADE simulation, the associated convergence of me- 
dian diameters is more pronounced there. 

In the coarse mode size range the difference in the 24 h 
mass size distributions is due to the inclusion of the HC1/C1 
equilibrium in MADE3. As we initialise the gas phase with- 
out HC1 (see Sect. 3), equilibration requires that some of the 
Cl initially evaporates from the particles (see Fig. 4). This 
reduction in Cl is responsible for the coarse mode particles’ 
mass loss in MADE3 with respect to MADE. 

4.1.2 Composition 

The temporal evolution of total aerosol species mass con- 
centrations in MADE3 and MADE is plotted in Fig. 4. 
In both MADE3 and MADE the tracers for mineral dust 
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Composition: MADE3 vs. MADE 



Figure 4. Temporal evolution of total aerosol species mass concen- 
trations in MADE3 (solid lines) and MADE (dotted lines). Abbre- 
viations are as follows: BC for black carbon and DU for mineral 
dust. Note that the particulate organic matter (POM) concentration 
is below 0.001 pgm -3 and therefore does not appear in the plot. 


(DU), black carbon (BC), particulate organic matter (POM, 
below 0.001 pgm -3 in our test case), and Na (as part of 
the sea spray tracer in MADE) always remain in the con- 
densed phase. Since we neglect particle mass sinks other than 
evaporation here, the results for these tracers are therefore 
identical. SO4 formation is faster in MADE 3 than in MADE 
due to the inclusion of H2SO4 condensation on the coarse 
mode particles. The most significant differences are seen in 
the H2O, Cl, NO3, and NH4 evolutions, where the Cl devi- 
ation was already described in the previous section and the 
loss of H2O in MADE 3 is due to the loss of Cl. 

The NH4 uptake in MADE 3 in the beginning of the simu- 
lation is coupled to SO4 uptake into the soluble Aitken mode 
particles. This process also occurs in the soluble accumula- 
tion mode, but only after NaNC>3 has been completely dis- 
placed by Na2SC>4 (~ 600 min; note that no Na is present 
in the Aitken modes in our test case). In EQSAM sodium 
ions and sulfate are neutralised first. NO3 therefore evapo- 
rates from the soluble accumulation mode particles because 
of the condensation of H2SO4 and subsequent replacement 
of NaNC>3 by Na2SC>4. When sulfate can no longer be neu- 
tralised by Na2SC>4 formation alone it becomes available for 
neutralisation by ammonium, leading to uptake of the lat- 
ter into the particles. This transition is visible as the kink 
in the MADE 3 NH4 curve (~ 600 min). The same applies 
to the MADE Aitken and accumulation modes (see kink at 
~ 300 min) but proceeds faster because coarse mode parti- 
cles are not a sink for the semi-volatile gases in MADE. 
This missing sink is also the reason for the second kink 
in the MADE NH4 evolution. As all the H2SO4 condenses 
on the fine particles in MADE, they eventually enter the 
sulfate rich regime (~ 1050 min). From this point on, sul- 
fate ions can bind less ammonium ions because EQSAM 
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Size distributions: MADE3 vs. PartMC-MOSAIC 
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Figure 5. Initial (grey, black) and final (light red, red) size distri- 
butions in MADE3 and PartMC-MOSAIC vs. dry diameter (i.e. ne- 
glecting aerosol water). Light colours (grey, light red) are used for 
number distributions (left vertical axis), full colours (black, red) for 
dry mass distributions (right vertical axis). MADE3 output is shown 
as solid lines, PartMC-MOSAIC output as crosses. The inset mag- 
nifies the coarse mode dry mass size distributions. 


then assumes that sulfate exists in the forms of (NH4)2S04, 
(NH4)3H(S04)2, and NH4HSO4 in the aerosol, whereas only 
(NH4)2S04 is considered in the sulfate neutral and sulfate 
poor regimes. 

NO3 is taken up by the coarse mode sea spray particles 
in MADE 3 via NaNC>3 formation despite the loss from the 
soluble accumulation mode. MADE, however, cannot rep- 
resent this process because it considers coarse mode par- 
ticles as passive and only a small amount of sea spray is 
present in the accumulation mode size range in our simu- 
lations. Hence, we see a continued increase in aerosol NO3 
content for MADE 3 , whereas this component is completely 
removed from the MADE aerosol. 

4.2 Comparison of MADE3 and PartMC-MOSAIC 

4.2.1 Size distributions 

In analogy to Sect. 4 . 1 . 1 , we first compare the size distribu- 
tions calculated by MADE 3 and PartMC-MOSAIC (Fig. 5 ). 
Note that the PartMC-MOSAIC results are averaged over 10 
runs, but the variability is less than the size of the crosses 
in the figure. Only at the large-diameter and small-diameter 
limits of the size distributions as shown here, is the variabil- 
ity higher because of the very few available computational 
particles. 

In the coarse mode size range (> 2 pm) the 24 h distri- 
butions of the two models agree very well. The disagree- 
ment between MADE 3 and PartMC-MOSAIC in the Aitken 
mode size range exposes a weakness of the modal approach 
with fixed mode widths. The Aitken mode becomes very 
narrow over the course of the PartMC-MOSAIC simulation. 


www.geosci-model-dev.net/7/1137/2014/ 


Geosci. Model Dev., 7, 1137-1157, 2014 


1150 

However, such narrowing cannot be simulated by MADE3, 
as the widths of its modes remain constant. The fast growth 
of the smallest particles (and slower growth of the larger 
Aitken mode particles) by condensation, for instance, can 
thus not be captured as accurately in MADE3. The 24 h 
MADE3 size distribution therefore contains more particles 
of very small diameters than the corresponding PartMC- 
MOSAIC distribution. 

Furthermore, although total number and mass concentra- 
tions of freshly emitted BC particles are the same in both 
models, the size distributions upon emission are different. 
PartMC-MOSAIC can use the original distributions from the 
emissions data set as given in Table 4 (see the “shoulder” to 
the right of the sharp peak in the 24 h number size distribution 
in Fig. 5). In MADE3, the particles are added to the wider 
modes, so that their contribution to the number size distribu- 
tion cannot be distinguished in Fig. 5. In addition, particle 
aging contributes to the less pronounced Aitken mode peak 
in MADE3. When the emitted BC particles acquire a coat- 
ing that surpasses the mass threshold of 10 %, they are trans- 
ferred to the mixed Aitken and accumulation modes. This 
leads to an increase of the median diameter of the mixed 
Aitken mode and to a reduction of the median diameter of 
the mixed accumulation mode. Hence, the two modes are no 
longer separately visible in the total number size distribution. 

Quicker growth of soluble Aitken mode particles in terms 
of both mass and diameter further adds to the shift of the peak 
in the MADE3 number size distribution with respect to the 
PartMC-MOSAIC distribution. That growth is due to water 
uptake on these particles that is predicted by EQSAM, but not 
by MOSAIC. Hence, some of the particles are renamed from 
the soluble Aitken to the soluble accumulation mode, so that 
the median (dry) diameter of the soluble accumulation mode 
is reduced. The latter mode thus also contributes to the wide 
peak in the MADE3 number size distribution. 

In conclusion, we see potentially significant differences 
between MADE3 and PartMC-MOSAIC in the size ranges 
of the fine particles. Such deviations had to be expected 
due to the simplifications and restrictions that come with 
the modal approach to represent particle size distributions. 
Despite these differences, simulation results with both mod- 
els agree well in the coarse mode size range (see also the 
next section and the size -resolved composition plots in Ap- 
pendix A). We are therefore confident that the coarse mode 
particle interactions that were newly introduced in MADE3 
are properly represented. 

4.2.2 Composition 

The evolution of total species mass concentrations gener- 
ally agrees well between MADE3 and PartMC-MOSAIC 
(Fig. 6 ). It should be noted that this agreement was achieved 
only after inclusion of the HC1 / Cl equilibrium in EQSAM, 
which leads to the decrease in Cl concentration and to the 
associated reduction of aerosol water content as described 
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Figure 6. Temporal evolution of total aerosol species mass con- 
centrations in MADE3 (solid lines) and PartMC-MOSAIC (dotted 
lines). Abbreviations are as follows: BC for black carbon and DU 
for mineral dust. Note that the particulate organic matter (POM) 
concentration is below 0.001 pg m - -’ and therefore does not appear 
in the plot. 


in Sects. 4.1.1 and 4.1.2. The difference in the mineral dust 
(DU) concentration is due to the very low number concen- 
tration of dust containing particles. The stochastic nature of 
PartMC-MOSAIC thus leads to a relatively large spread of 
the DU concentrations in the ten simulations (not shown). 
Since the MADE3 result falls within the range of simulated 
values this deviation does not impair the overall agreement. 

NO 3 is taken up more quickly in MADE3 than in PartMC- 
MOSAIC due to the assumption that equilibrium is attained 
during each time step. Note that the flux limit described in 
Sect. 2. 1 .2 under “Gas-particle partitioning” is never reached 
in our test case. In addition, there is more NO 3 partitioning 
to the condensed phase in MADE3 as it displaces Cl from 
the particles, leading to a slightly lower Cl content of the 
MADE3 aerosol. These differences can be explained by the 
different chemistry codes. While MOSAIC allows for the co- 
existence of NaCl and NaN 03 at arbitrary Na concentrations, 
NaCl can exist in EQSAM only when all the available nitrate 
has been bound to the sodium ions. 

SO 4 uptake is slightly slower in PartMC-MOSAIC than 
in MADE3. As H 2 SO 4 condensation is limited by gas phase 
diffusion, the difference is due to the different assumptions 
for the accommodation coefficient: cv|[ 2 so 4 = 1 in MADE3, 
SO 4 =0.1 in PartMC-MOSAIC. Hence, the H 2 SO 4 flux 
to the particles is greater in MADE3 than in PartMC- 
MOSAIC (see Eq. 9). This was confirmed by a sensitivity 
simulation, in which we set aj^scq =0.1 in MADE3. 

The only qualitative difference between the two models in 
terms of composition evolution is in the NH 4 concentration. 
It is due to the different treatments of activity coefficients 
by EQSAM and MOSAIC (for details, see Metzger et al., 
2002; Zaveri et al., 2005b). While the activity coefficient in 
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MOSAIC allows some condensed phase NH 4 CI to be pro- 
duced in the accumulation mode size range, EQSAM pre- 
dicts that this compound will not form in our test case. 

In MADE3 ammonium partitions only to the fine particles 
because there are not enough anions available in the coarse 
modes to neutralise it. In the beginning of the simulation 
NH 4 is formed only in the smallest particles as (NELj^SCU. 
These particles eventually grow into the accumulation mode, 
where NO 3 is taken up as NaN 03 because HNO 3 is available 
more quickly than H 2 SO 4 . However, SO 4 from H 2 SO 4 con- 
densation eventually displaces the NO 3 in the MADE3 ac- 
cumulation mode particles. When all NO 3 has left these par- 
ticles (~ 600 min) the additionally condensing SO 4 is avail- 
able for (NH 4 ) 2 SC >4 formation. Thus, the NH 4 concentration 
rises more quickly thereafter. 

5 Summary and conclusions 

We have presented MADE3, a modal aerosol submodel with 
nine lognormal modes for use within MESSy as part of the 
AC-GCM EMAC. The modes represent the three classes of 
fully soluble, insoluble, and mixed particles in each of three 
size ranges, namely the Aitken, accumulation, and coarse 
mode size ranges. The submodel has been developed on the 
basis of its predecessors MADE (Ackermann et al., 1998; 
Lauer et al., 2005) and MADE-in (Aquila et al., 2011). It was 
extended by inclusion of coarse mode particle interactions 
with condensable trace gases and with other particles, and 
by inclusion of the gas-particle partitioning of chlorine. We 
compared the new submodel in a box model setup to its pre- 
decessor MADE and to the state-of-the-art particle-resolved 
aerosol box model PartMC-MOSAIC (Riemer et al., 2009; 
Zaveri et al., 2008), which we used as a reference here. For 
the comparison we designed and discussed an example test 
case representative of clean marine boundary layer condi- 
tions with added shipping emissions. This setup was chosen 
because coarse mode particle interactions potentially play an 
important role in such an environment. 

We obtained similar results with MADE3 and MADE, but 
there were differences especially in the coarse mode size 
range. We expect to find less Cl and more NO 3 in the aerosol 
particles on global average when switching from MADE to 
MADE3. The evaporation of some of the Cl to the gas phase 
(as HCI) may entail differences in aerosol water content if 
the lost particulate Cl is not fully replaced by NO 3 . This, in 
turn, may affect the prediction of cloud condensation nuclei 
(CCN) concentrations and aerosol optical depth calculations. 
Total aerosol mass concentrations may also differ whenever 
sea spray particles dominate the aerosol mass concentration. 


Comparing MADE3 to PartMC-MOSAIC, we found good 
agreement in tenns of total aerosol composition evolution 
and coarse particle (> 2 pm) size distribution predictions, de- 
spite some potentially significant differences in the size dis- 
tributions of fine particles. It is important to note that the par- 
ticle size distribution is one of the main factors that govern 
the conversion of aerosol particles to cloud droplets. In our 
test case, MADE3 results show a particularly large deviation 
from the PartMC-MOSAIC results in the size range where 
this activation primarily takes place. Only the 3-D model 
application of MADE3, i.e. the inclusion of the processes 
omitted here and the consideration of many different combi- 
nations of environmental parameters and aerosol population 
states, will allow us to check against observations whether 
this deviation can be generalised and whether it could lead 
to a systematic bias in aerosol-cloud interactions. According 
to the results of a comparison of a modal aerosol submodel 
with a sectional one within the same global chemical trans- 
port model (Mann et al., 2012), CCN concentrations may be 
overestimated by the modal approach in certain regions. 

Nevertheless, since MADE has been extensively and suc- 
cessfully evaluated as part of different AC-GCMs and chem- 
istry transport models (e.g. Lauer et al., 2005; Ochoa et al., 
2012; Zhao et al., 2013; Righi et al., 2013), we are confident 
that MADE3 is also suitable for use with EMAC. Consider- 
ing the similar results obtained with MADE3 and MADE in 
the box model test case that was drawn from an actual 3-D 
model run, we expect similar performance for MADE3 in 
the 3-D model as well. From the comparison with PartMC- 
MOSAIC we conclude that improvements in the representa- 
tion of coarse mode aerosol particles and total aerosol com- 
position are likely when switching from MADE to MADE3. 
A corresponding evaluation within the MESSy framework 
by means of comparison with observational data will be the 
subject of a follow-up study. 
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Appendix A: Evolution of size-resolved 
aerosol composition 

In Fig. A1 we show the initial and final states of the aerosol 
population as simulated by MADE3 and PartMC-MOSAIC. 
The figure illustrates the evolution of the aerosol composi- 
tion under additional consideration of the size distribution. 
Composition is shown in a size-binned representation to fa- 
cilitate the comparison of the individual panels. The top left 
panel contains the same data as Fig. 2 (note that the mass 
fractions of mineral dust are so small that they are not visible 
in Fig. Al). 

The lower row plots in Fig. Al, representing the aerosol 
state after 24 h of simulated time, show the same general fea- 
tures: higher SO 4 fractions in the fine particles, predominant 
NO 3 partitioning to large rather than small particles, and 
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notable BC concentrations only in fine particles. Despite this 
agreement, one can also see that the modal approach leads to 
distribution of the components over wider size ranges. Since, 
due to their width, the modes contain particles of a broad 
range of sizes, this is inevitable in modal models. 

The smallest particles do not take up water in PartMC- 
MOSAIC because they are assumed to be dry initially (Za- 
veri et al., 2005a) and the deliquescence relative humidity 
of (NH 4 ) 2 S 04 is higher than the environmental relative hu- 
midity specified in our experiment (0.771). Conversely, in 
EQSAM, these particles do take up water due to the pres- 
ence of small quantities of other components that reduce the 
particles’ deliquescence relative humidity. 
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MADE3 wet mass size distribution (t = 0) 



MADE3 wet mass size distribution (t = 24h) 
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PartMC-MOSAIC wet mass size distribution (t = 0) 
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Figure Al. Aerosol mass size distributions and size-resolved composition at t = 0 (top panels) and t = 24 h (bottom panels) for MADE3 
(left panels) and PartMC-MOSAIC (right panels). Total mass size distributions are shown in the MADE3 panels as black lines (calculated as 
the sums over the modes), in the PartMC-MOSAIC panels as black crosses (left vertical axes). The coloured bars show the contributions of 
the individual species (POM = particulate organic matter, BC = black carbon) to the respective mass concentrations (right vertical axes). In 
the PartMC-MOSAIC plots particles were binned into size sections, and the bars show the average compositions in these sections. For the 
MADE3 plots the binned mass concentrations were computed as the sums (across all modes) over the fractions of particles that fell within 
the same size sections. Note that while the right vertical axes are logarithmic, the species fractions in the bars add up linearly to the total 
mass concentrations (i.e. the axes only apply to the total mass concentrations in each bin, but not to the individual contributions). 
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Code availability 

The Modular Earth Submodel System (MESSy) is continu- 
ously further developed and applied by a consortium of insti- 
tutions. The usage of MESSy, including MADE3, and access 
to the source code is licensed to all affiliates of institutions 
which are members of the MESSy Consortium. Institutions 
can become members of the MESSy Consortium by sign- 
ing the MESSy Memorandum of Understanding. More in- 
formation can be found on the MESSy Consortium Website 
(http://www.messy-interface.org). 
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